Structural properties of screened Coulomb balls 



in 
o 
o 

(N 

< 
o 

m 

Oh: 



o ■ 

• I— I . 



> 

(N 

^ ' 

00 

o 

^ ; 
o : 



X 



M. Bonitzi, D. Block^B O. Arp^, V. Golubnychiy^ , H. Baumgartner^, P. Ludwig^, A. Piel^, and A. Filinov^ 

^ITAP, ^lEAP, ChnsUan-Albrechts-Umversitat zu Kiel, D-24098 Kiel, Germany 

(Dated: February 2, 2008) 

Small three-dimensional strongly coupled charged particles in a spherical confinement potential 
arrange themselves in a nested shell structure. By means of experiments, computer simulations and 
theoretical analysis, it is shown that their structural properties depend on the type of interparticle 
forces. Using an isotropic Yukawa interaction, quantitative agreement for shell radii and occupation 
is obtained. 
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The recently discovered Coulomb balls are an inter- 
esting new object for studying strongly coupled systems. 
Coulomb balls consist of hundreds of micrometer sized 
plastic spheres embedded in a gas plasma. The plastic 
spheres attain a high electric charge Q of the order of 
several thousand elementary charges and arrange them- 
selves into a highly ordered set of nested spherical shells 
with hexagonal order inside the shells. Coulomb balls 
are a special form of 3D-plasma crystals ,3„ . The 
formation of ordered clusters with nested shells was also 
observed in laser-cooled trapped ion systems, e.g. [^|^, 
and is expected to occur in expanding neutral plasmas 

m 

The same kind of ordering was found in molecular dy- 
namics (MD) simulations, e.g. 9, 10, 11] and references 
therein. In particular, the transition to the macroscopic 
limit [T3 | , the symmetry properties of the individual 
shells including a Voronoi analysis ^3 and metastable in- 
trashell configurations [Til fl^ have been analyzed. Very 
large systems of trapped ions show a transition to the 
crystal structure of bulk material, which was measured 
by laser scattering Il5|. 

Although the shell structure of ion crystals is quite 
well understood in terms of simulation results, these sys- 
tems do not yet allow for systematic experimental studies 
of the structure inside the shells and the detailed occu- 
pation numbers of individual shells. The advantage of 
studying Coulomb balls is the immediate access to the 
full three-dimensional structure of the nested shell sys- 
tem by means of video microscopy. The tracing of each 
individual particle is ensured by the high optical trans- 
parency of the system, which results from particle diam- 
eters of typically 5 /im at interparticle spacings of 500 
/im. Compared to atomic particles, the very high mass 
of the microparticles used here slows down all dynamic 
phenomena to time scales ranging from 10 ms to sec- 
onds. Therefore, studies of Coulomb balls complement 
investigations of ion crystals, where dynamical studies 
are difficult. 

Coulomb balls are in a strongly coupled state, i.e. 
the Coulomb coupling parameter, T — ^ , where a 
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is the mean interparticle distance, attains large values 
(r 3> 100). Contrary to ion crystals, where the parti- 
cles interact via the pure Coulomb force, the micropar- 
ticles in a Coulomb ball are expected to interact by a 
Yukawa type potential, Vdd = Q e~'"'^°/r, which was ver- 
ified experim<;ntally in complex plasmas M 

Therefore, 

Coulomb balls are characterized by two parameters, the 
coupling parameter T and the Debye shielding length of 
the plasma A^i. It is the intention of this paper to study 
the influence of shielding on the structure of Coulomb 
balls, in particular, to pin down the differences from sys- 
tems with pure Coulomb interaction. This will be done 
by comparing computer simulations with experimental 
results. At the same time, a study of spherical crys- 
tals with Yukawa interaction opens up an interesting new 
field which in a natural way bridges the gap between the 
above mentioned theoretical investigations of finite size 
Coulomb systems and the theory of macroscopic Yukawa 
plasmas, e.g. [r^lH|. 

Experiment. The experiment is described in detail in 
Refs. so here we only summarize the main re- 

sults from a systematic investigation of 43 Coulomb balls 
consisting of 100 to 500 particles. All Coulomb balls were 
trapped under identical experimental conditions. All of 
them had a spherical shape and their diameter was in 
the range of 4-5 mm. A typical experimental result for 
a cluster and its shell structure is shown in the left part 
of Fig. n In all 43 Coulomb balls a similar shell struc- 
ture was observed and the shell radii Rg and the shell 
occupation numbers Ns were measured. Further, from 
the pair correlation function the typical mean interpar- 
ticle distance was derived, which for all A'' was close to 
a ~ 0.6 mm. The mean intershell distance d was found 
to be c? = (0.86 ± 0.06)a, which yields a ratio d/a ~ 0.82 
in good agreement with local icosahedral ordering |^. 
An important experimental result is that the intershell 
distance is constant over the whole Coulomb ball and 
implies a constant average density of particles and ions, 
which is equivalent to a parabolic electric potential well 
used for the simulations below. A more detailed experi- 
mental verification of the parabolic confinement well will 
be described elsewhere [l^l- A different case with 'self- 
confinement' of a dust cloud in a strongly anharmonic 
potential was recently discussed in [20l |. 

Simulations. For a theoretical explanation of the ex- 
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FIG. 1: Radial particle distribution for A'^ = 190 given in 
cylindrical coordinates. Left: experiment P], right two figures: 
simulation results with Coulomb (k = 0), and Yukawa (k = 1) 
potential. The length unit in the right two figures is Toc, given 
byEq. 

perimental results we have performed molecular dynam- 
ics (MD) and thermodynamic Monte Carlo (MC) simu- 
lations using the Hamiltonian 

We assume that the Coulomb balls consist of particles 
with same mass and charge and that a stationary state 
is reached close to thermodynamic equilibrium. Further- 
more, the observed isotropic particle configuration sug- 
gests to use an isotropic interaction potential. Screen- 
ing effects are included in static approximation using 
Debye(Yukawa)-type pair potentials Vdd given above. In 
the simulations we use dimensionless parameters, with 
lengths given in units of the ground state distance of 
two particles, tqc, defined in Eq. ||2Jl, hence in this paper 

= fQc/^D- In experimental papers, k = a/Xo is often 
used. In accordance with the experiment on Coulomb 
balls and previous experiments and simulations on 
ion crystals 17|, we use a screening- independent confine- 
ment potential Uc{r) — muj^ ■ As a result, in our 
model the configuration of the Coulomb balls is deter- 
mined by three parameters: particle number N , screen- 
ing parameter n and temperature T. Since experimen- 
tal plasma densities and temperatures are not precisely 
known, we have performed a series of calculations for dif- 
ferent values of k and T. Furthermore, a wide range of 
particle numbers, up to A'^ = 503, has been analyzed. 

Results. Consider first the theoretical ground state 
configurations (T = 0) in the case of Coulomb interac- 
tion, K = 0, which were obtained by classical MD simula- 
tions using an optimized simulated annealing technique 
[Tl| . Using about 1000 independent runs for each value of 
N ensured that the ground state is reached. In addition, 
we have performed MC simulations in the canonic ensem- 
ble with a standard Metropolis algorithm, which allows 
for a rigorous account of finite temperature effects. Both 
simulations yield identical configurations at low temper- 
ature. Fig. n shows a comparison of MD simulation and 
experiment for the case of iV = 190 particles. In both 



cases four concentric spherical shells are observed, which 
are the result of a balance between confinement potential 
Uc and interparticle repulsion Vdd- 



4 



3 







— K-0 

-- K-0.6 
























O 












i 
















< 




A 




1 









FIG. 2: Experimental (symbols) and simulation (lines) results 
for the shell radii of three-dimensional Coulomb balls in units 
of the mean interparticle distance. 

For a more detailed quantitative comparison between 
experiment and simulation we analyze the dependence of 
the shell radii Rs on the cluster size iV (Fig. [^J. The 
interparticle distance o serves as a common length scale 
as it is accessible in experiment and simulation. There 
is an overall increase cx N^l'^ of the experimental Rs for 
all shells and all 43 analyzed clusters. Exceptions occur 
around values of N where new shells emerge. The same 
behaviour is obtained from the MD simulations. Without 
any free parameter a very good agreement of experimen- 
tal radii and Coulomb MD results (full lines) is observed, 
in particular concerning the absolute values, the slope 
and the equidistance of the shells. Further, this results 
holds even in case of a Yukawa potential if k is small 
(dashed lines in Fig. |2J|. These findings already imply 
two things: First, the approach to use an screening in- 
dependent parabolic confinement potential to model the 
experiments is justified. This is a marked difference from 
a situation with self-confinement, which does not yield 
equidistant shells jl^]. Second, the general scaling of the 
shell radii of weakly shielded Coulomb balls cx N^l'^ is the 
same as for pure Coulomb systems, such as ion crystals. 

However, a marked difference between experiment and 
simulations of pure Coulomb systems is observed for the 
shell population numbers Ni . . . N4. Table ^ shows the 
shell population numbers for various screening param- 
eter K of a Coulomb ball with N = 190 as obtained 
from MD simulations and experiment. Clearly, the MD 
results yield systematically more particles in the outer 
part of the cluster for k = than observed in exper- 
iment. Further, Tab. ^ shows that, with increasing k, 
particles move from the outer shell inward. Interestingly, 
for K = 0.58... 0.63, the simulations yield exactly the 
same shell configuration as the experiment. Hence, the 
different population numbers may be attributed to the 
influence of screening and hence weakening of the inter- 
action potential. 

To investigate this in more detail, this comparison was 
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extended to all 43 Coulomb balls. Due to their dif- 
ferent size and even different number of shells the sys- 
tematic differences in shell population of Coulomb and 
Yukawa systems can be studied comparing the experi- 
mental results and systematic MD-simulations. The re- 
sult is shown in Fig. 13 For all systems Coulomb, Yukawa 
as well as the experimental data the shell population of 
all shells shows an almost linear behavior as a function of 
iV^/^. However, the experimentally obtained shell pop- 
ulation of the outermost shell N4 is always smaller than 
those of a Coulomb system (solid line) and the inner 
shells show a systematically higher population. From the 
Yukawa MD-simulations (dashed lines) this tendency is 
confirmed. It is clearly found that particles move to in- 
ner shells with increasing k. Hence, the finding discussed 
for the Coulomb balls with N = 190 in Tab. |l| holds 
generally. This tendency reflects the fact that, from an 
energetic point of view, the higher population of inner 
shell becomes less costly due the shielding than the oc- 
cupation of the outermost shell, where the confinement 
by the trap dominates the potential energy. 




FIG. 3: Experimental (symbols) and simulation (lines) re- 
sults for the shell population of three-dimensional Coulomb 
clusters. 

In more detail, we find that the outermost shell ex- 
hibits the largest absolute change with n and, there- 
fore, it is best suited for a detailed comparison with 
the experimental data, see Fig. |21 From a best fit to 
the experimental data, we find a screening parameter 
^EXP _ Q g2 ± 0.23. An independent analysis for the 
other shells confirms this result, e.g. the third shells, 
yield k^-^^ = 0.58 ± 0.43. Determining the mean inter- 
particle distance a from the first peak of the pair distri- 
bution function k^^^ translates into an average Debye 
length Xo/a = 1.54 ± 0.7. Furthermore, as one can see 
in the right hand part of Fig. an increase of k leads 
to compression of the entire cluster, which is due to the 
reduction of the potential Vdd- The fact that more and 
more particles move from the outer shells inward has the 
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TABLE I: Experimental (last column) and theoretical shell 
configuration of the Coulomb ball N = 190. Ni . . . N4 denote 
the particle numbers on the i-th shell beginning in the center. 



consequence that closed shell configurations are already 
reached for a smaller number N* of total particles com- 
pared to N* in the Coulomb case. While for k = 0, the 
first closed shell is found at 12 particles, for k > 4.7 the 
ground state of a cluster with 12 (and 11 as well) par- 
ticles contains one particle in the center and N* = 10. 
For K = 0.6 closure of the 2nd to 4th shell is observed 
for N* = 54, AT* = 135, = 271_,^compared to the 
Coulomb case where N^^ = 57,60 [13, N^c = 154 [l| 
and N^^ = 310, cf. also Fig.H 

After analyzing the shell populations we now consider 
the shell width. The larger roughness of the shells in the 
experiments, cf. Fig.^ is attributed to small anisotropics 
of the experimental confinement and finite depth resolu- 
tion of the imaging equipment as well as temperature 
effects. While the measurements are at room tempera- 
ture, the simulations refer to T = 0. Therefore we have 
analyzed the influence of temperature on the shell radii 
and populations with the help of MC simulations. From 
the results we conclude that the effect of temperature on 
the shell configurations Ns is negligible for k = 0.6. At 
constant finite T we find that an increase of k leads to 
a reduction of shell roughness. Contrary to that, a tem- 
perature increase at elsewhere fixed parameters in fact 
leads to a roughening of the shells proportional to VT 
for the outer shell and an even stronger effect for the in- 
ner shells. This tendency will become evident from the 
analytical results below. 

Analytical results. The main influence of screening on 
Coulomb balls is readily understood with the help of an- 
alytical results, which can be found for N = 2. First, the 
ground state distance ro(«;) follows from minimizing the 
potential energy U in Eq. 



pK.ro „3 
6 '0 

1 -I- Kro 



— 'Oc- 



(2) 



Eq. ^ yields the two-particle distance, Tqc, in an un- 
screened system as a function of rg and is easily inverted 
numerically '21] . The ratio tq /tqc is always smaller than 
unity and monotonically decreasing with k, thereby con- 
firming the above observation of screening-induced com- 
pression of the Coulomb balls. Second, we analyze the 
cluster stability by expanding the potential U in terms of 
small fluctuations, y = r — vq, around the ground state, 
up to second order: U{r)-U{ro) = ^U"{ro)y^ = ^fl^y^. 
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This defines an effective local trap frequency 

= 30.2 f 1 + ^T^) = -^f2{n), (3) 
\ 3 1 + Kro / m Tq 

/2H = e-«'-° (1 + Kro + ^2^2/3) , 

which allows us to estimate the width of the Coulomb 
ball shells. Finally, we compute the variance of the par- 
ticle distance fluctuations, ct^, for particles in a parabolic 
potential with frequency Q at temperature T and ob- 
tain = AksT / [m^^). This allows for two interesting 
conclusions, which explain the above MC results: At con- 
stant screening, the shell width grows with temperature 
as "v/T while screening reduces the shell width. One might 
be tempted to conclude that increased screening makes 
particle transitions between neighboring shells less likely 
and thus stabilizes the cluster against melting. However, 
the opposite is true, because screening also reduces the 
distance between shells which is of the order of rp. The 
relative importance of both tendencies can be discussed 
in terms of the relative distance fluctuations, a critical 
value of which determines the onset of radial melting 
(Lindemann criterion). 

= r* = r2/2M . (4) 

'0 ^ 2 

Ur is related to an effective coupling parameter, which 
depends on the interaction strength of two trapped par- 
ticles - via the Coulomb-type coupling parameter, F2 = 
/ (kBTro), and on the screening strength - via the 
function f2{n)- /2 monotonically decreases with k (ur in- 
creases), thus screening always destabilizes the Coulomb 
balls. 

Finally, these analytical results are closely related to 
those for macroscopic homogeneous Yukawa systems, e.g. 
[TtI ITsII. This limit is recovered by replacing, in 
by the mean interparticle distance a at a given density n, 



a = (3n/47r)"'^/^. Then the local trap frequency becomes 
fi^ —> ujp^f2{K), showing that, in a Coulomb system, Q 
approaches the dust plasma frequency LOpd whereas, in 
the case of screening, the result is modified by a factor 
f2{n) nil- Also, the effective coupling para meter Fj is 
in full analogy to the macroscopic result p!q . 

In summary, we have presented a combined experimen- 
tal, numerical and theoretical analysis of small spheri- 
cal charged particle clusters. The excellent experimen- 
tal accessibility of these systems has been demonstrated. 
The structure of these clusters deviates from models 
with pure Coulomb interaction and requires the inclu- 
sion of static screening. For the particle number range 
TV = 100 . . . 500, comparison with the MD and MC sim- 
ulations has allowed us to determine the screening pa- 
rameter averaged over the clusters asA^/assLS. These 
Coulomb balls are representative for finite Yukawa sys- 
tems, combining shell properties observed in spherical 
Coulomb clusters with screening effects found in Yukawa 
plasmas. Since the shell occupation numbers have now 
been critically analyzed, our results confirm earlier con- 
clusions about the shell structure of ion clusters, where 
such an analysis was not accessible yet. The results are 
relevant for other strongly correlated charged particle 
systems, such as crystal formation of droplets in expand- 
ing laser produced plasmas, where shielding becomes im- 
portant. 



Acknowledgments 

This work is supported by the Deutsche Forschungsge- 
meinschaft via SFB-TR 24 grants A3, A5 and A7 and, in 
part, by DLR under contract 50WM0039. We acknowl- 
edge discussions with W.D. Kraeft and M. KroU's assis- 
tance in conducting the experiments. 



[1] O. Arp, D. Block, A. Piel, and A. Melzer, Phys. Rev. 

Lett. 93, 165004 (2004) 
[2] J.B. Pieper, J. Goree, and R.A. Quinn, Phys. Rev. E 54, 

5636 (1996). 

[3] M. Zuzic et al. Phys. Rev. Lett. 85, 4064 (2000) 

[4] Y. Hayashi, Phys. Rev. Lett. 83, 4764 (1999) 

[5] D.J. Wineland, J.C. Bergquist, W.M. Itano, J.J. 

Bollinger, and C.H. Manney, Phys. Rev. Lett. 59, 2935 

(1987) 

[6] M. Drewsen et al., Phys. Rev. Lett. 81, 2878 (1998) 

[7] T. Pohl, T. Pattard, and J.M. Rost, Phys. Rev. Lett. 92, 

155003 (2004) 
[8] T. Killian, Nature 439, 815 (2004) 

[9] R.W. Basse, and V.V. Avilov, Phys. Rev. A 44, 4506 
(1991) 

[10] K. Tsuruta, and S. Ichimam, Phys. Rev. A 48, 1339 
(1993) 

[11] P. Ludwig, S. Kosse, and M. Bonitz, Phys. Rev. E 71, 



046403 (2005) 

[12] H. Totsuji et al., Phys. Rev. Lett. 88, 125002 (2002) 
[13] J.P. Schifler, Phys. Rev. Lett. 88, 205003 (2002) 
[14] O. Arp et al., J. Phys. Conf. Series 11, 234 (2005) 
[15] W.M. Itano et al.. Science 279, 686 (1998) 
[16] U. Konopka, G.E. MorfiU, and L. Ratke, Phys. Rev. Lett. 
84, 891 (2000) 

[17] D.H.E. Dubin, and T.M. O'Neill, Rev. Mod. Phys. 71, 
87 (1999) 

[18] V.E. Fortov, et al., Phys. Rev. Lett. 90, 245005 (2003) 

[19] O. Arp, D. Block, and A. Piel, to be published 

[20] H. Totsuji, C. Totsuji, T. Ogawa, and K. Tsuruta, Phys. 

Rev. E 71, 045401 (2005) 
[21] A useful analytical approximation for ro in Eq. ((2J as a 

function of roc is x = Xc + '°J2"^3^^''^3 " ^c{^ + Xc), where 

X — nro, which has an error of less than 1%, for Xc = 

KTOc < 1.5. 



[22] This result differs sfiglitly from tlie exact macroscopic tliis replacement in tfie two-particle expression 

result 13] [by the coefficient 1/3 instead of 1/2 in the [23] H. Baumgartner, W.D. Kraeft, and M. Bonitz, to be 
last term in /2] which is a consequence of performing fished 



